Decoherence of a qubit by non-Gaussian noise at an arbitrary working point 
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The decoherence of a qubit due to a classical non-Gaussian noise with correlation time longer than 
the decoherence time is discussed for arbitrary working points of the qubit. A method is developed 
that allows an exact formula for the phase memory functional in the presence of independent random 
telegraph noise sources to be derived. 



I. INTRODUCTION 

Reducing the decoherence induced by interaction with 
the environment is one of the major challenges in the 
practical implementation of quantum computing. In par- 
ticular, for solid state qubits this seems to be the most 
important issue. In this paper we discuss the effect of 
noise with a correlation time that is long compared to 
the decoherence time of the qubit. The work is directly 
motivated by experiments- 1,2 on Josephson charge qubits 
(JCQ), but due to their general nature they can also be 
relevant to many other systems. The JCQ is built around 
a small superconducting grain connected to a supercon- 
ducting reservoir by a Josephson junction. By means of 
a capacitively coupled gate voltage one can control the 
number of Cooper pairs on the small grain. Because of 
the Coulomb blockade there will be a preferred number of 
Cooper pairs except at special degeneracy points where 
the energies of states with n and n + I coincide. The de- 
generacy is lifted by the Josephson coupling, giving the 
usual level anticrossing picture sketched in Fig. 2] If the 
gate voltage is never too far from the degeneracy point 
we can ignore transitions to other levels and the device 
can be regarded as a two level system or qubit. 

Consider noise in this system originating from fluctua- 
tions in the electric potential of the grain. This could be 
either due to fluctuations in the gate voltage source, or 
to fluctuating charges in the environment. Both would 
correspond to fluctuations of the working point position 
along the horizontal axis of the figure, with a correspond- 
ing change in the energy. The experiments^ where con- 
ducted with the gate voltage away from the degeneracy 
point, where to a good approximation the energy is a lin- 
ear function of the potential. This case is covered by the 
theorji^ for the case of free induction and narrow distri- 
bution of coupling constants and for both free induction 
and echo at arbitrary distribution of coupling constants, 
respectively. 

Using a modified circuit Devoret et. alSk were able to 
work at the degeneracy point where to first order there 
is no change in energy when the potential is changed. 
The idea is that this would make the device less sensitive 
to electrostatic noise. This is then called the optimal 
working point. At this point the expectation value of 



the charge is the same for both states. The purpose of 
this paper is to extend the theory2i4 to this situation and 
to see how dephasing is changed as the optimal point is 
approached. 

The paper is organized as follows. In section [H] the 
model for the qubit interacting with the noise is pre- 
sented. In section lIlTl we discuss a short time expansion 
that illustrates in a simple way the interplay of the effect 
of several uncorrelated fluctuators at the optimal point. 
The main purpose of this section is providing a qualita- 
tive understanding, as the results are contained within 
the full solution presented in section ITVl 



II. MODEL 

Consider the Hamiltonian of the Josephson qubit, c. f. 
withRef.H 



H =^Acr z - ^Ejcr x 



(1) 



g \ y g - V° pt ) /e, where E c is the charging 



Here A = E c C g (V & 
energy, (V g — V° pt ) is the deviation in gate voltage from 
the optimal working point, C g is the gate capacitance and 
e the electron charge. The energy levels as function of 
the gate potential are shown in Fig.^ In principle, both 
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FIG. 1: Energy levels of Josephson qubit as function of the 
gate potential. 

A and Ej can fluctuate in time, but for clarity we will 
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only consider the noise in A. Noise in Ej was considered 
in Ref. 6 and can be treated in a similar way. 

We will describe fluctuations in A by a random time- 
dependent quantity v(t). The Hamiltonian of the qubit 
is then 



v(t)]a z - ^Eja x 



(2) 



One can express v(t) through fluctuations of the effec- 
tive gate voltage, 5V g (t) as v(t) = E c C g 8V g (t)/e. One 
of the sources of such fluctuations is hopping background 
charges i&^i 7 . We briefly recall the model of Refs. It 
is assumed that there are dynamic defects characterized 
by two metastable states. They can be, for example, 
traps near the electrodes able to capture electrons from 
the electrodes and then re-emit them, or pairs of traps 
with one electron bouncing between them. Each dynamic 
defect is represented as a classical fluctuator producing 
random telegraph noise. That is, the state of each fluc- 
tuator is represented by a random function \(t) which 
switches between the two values, ±1/2. The probability 
to switch n times between theses states during the time 
t is assumed to be given by the Poisson distribution 



Pn{t) 



{itT - 



It 



(3) 



Here 7 is the characteristic switching frequency of the 
given fluctuator. This means that the switching events 
occur independently of each other. The fluctuation v(t) 
is a sum of the contributions of different fluctuators, 



(4) 



Different fluctuators are assumed to be statistically inde- 
pendent, 
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-2-Vi|ti 



(5) 



This is actually the simplest model since it is assumed 
that the probabilities to jump between the states in two 
directions are the same. As a result, equilibrium popu- 
lations of the states are equal. This assumption can be 
relaxed- without change of the physical conclusions. In- 
deed, since only fluctuators with energy splittings < kT 
contribute to the decoherence, the equilibrium popula- 
tions of the efficient states do not essentially differ. 

Several paper s 3 ' 4 ' 7 have addressed the dephasing in the 
presence of a large number of fluctuators with a broad 
range of switching rates, 7. This situation is relevant for 
systems showing 1// noise. We believe, however, that the 
experiments with "good" qubits, e. g., Ref. Q], are best 
understood in terms of a small number of fluctuators. 
Also, if one is to achieve quantum computation this is 
necessary as a large number of fluctuators most probably 
will lead to large decoherence and ruin the device. We 
will therefore here focus on the effect of a small number 
of fluctuators. 



The noise affects the qubit in two ways: (i) it causes 
shifts in the energy levels of the two states and thereby 
introduces a random contribution to the relative phase of 
the two states (dephasing) and (ii) it causes transitions 
between the two states leading to energy relaxation. Let 
us for the moment concentrate on the first effect, de- 
phasing. The important quantity that we study is the 
phase-memory functional 



exp 



dt'(3(t')SE({ Xi (t')}) 



(6) 



where 5E({xi(t')}) is the shift in the energy splitting 
of the two levels caused by the fluctuators, while (3{t) 
is some function which depends on the qubit manipula- 
tion procedure, as will be explained below. The phase- 
memory functional describes the relative phase picked up 
during time evolution by one state of the qubit relative 
to the other. 

Diagonalizing Eq. we get the eigenenergies 



E± = ± 



(J) 



If the working point is far from the optimal one, A 3> Ej, 
we can neglect the Josephson energy, and coupling of the 
fluctuators to the qubit is linear, 



E± 



Ejll 



(8) 



An exact formula for the phase memory functional for a 
linear coupling in the case of a single fluctuator coupled 
to the qubit was derived in Refs. yEB In this limit, be- 
cause of the linearity E in v, the phase memory functional 
for any number of fluctuators is found by simply multi- 
plying the phase memory functionals of different fluctu- 
ators. In the regime of exponential decay this procedure 
corresponds to simply adding the decay rates and the 
resulting decay rate is represented by the average over 
the distributions of fluctuator parameters/ The optimal 
point, A = 0, was studied in various approximations and 
numerically in Refs. Ijjla. 

The aim of the present work is to derive an exact for- 
mula for the phase memory functional of a single fluctu- 
ator similar to that in Refs. 00 that is applicable at an 
arbitrary working point. Averaging over a large number 
of fluctuators is in the general case not as simple as in 
the linear regime since the phase memory functional is 
not the product of individual factors for each fluctuator. 
However, we will extend the analysis to a small num- 
ber of fluctuators, which we believe is relevant to qubit 
experiments^ 

To make the main idea clear let us first study the ex- 
pansion at short times, yt -C 1. Firstly, this expansion 
provides some insight; secondly, the short-time situation 
can be most important for realistic qubits since at long 
times the phase memory functional has already decayed 
to a very low value. 
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III. PHASE MEMORY FUNCTIONAL AT 
SMALL TIMES, jt < 1 

We expand the energy in the limit 



to obtain 



v < E 



E± = ±-[ E 



A 2 



Av 



v 

2E* 



(9) 



where E* = Eq/Ej. When averaging the phases we sub- 
tract the initial values so that we only get what comes 
from fluctuator jumps. Thus \& = (e 1 ^') where 



1 



m = t / PVW 



E VjVj 
1e* 

ij 



[Xi(t')Xj(t') 



x°x°] 



(10) 



is the random phase shift. Here %® = Xi(0) while the 
average is calculated over random telegraph processes in 
the fluctuators. 

If there is only linear coupling, each fluctuator appears 
in only one term in the exponent. Since we assume the 
fluctuators to be statistically independent the average 
can be split into a product of averages over each indi- 
vidual fluctuator. With the quadratic term included this 
is no longer possible. Here we will study the cases of 
one and two fluctuators. This gives some insight to the 
general structure, and it is also the most relevant for real- 
izing a working qubit since a large number of fluctuators 
would destroy the working of the qubit. 

a. Free induction signal: Let us assume that (3{t) = 
9(t) where 6(t) is the Heaviside unit step function. This 
assumption corresponds to the free induction signal. 
With only one fluctuator the quadratic term is identi- 
cally zero, and we have 



A 



2Enh 



(11) 



Let us evaluate the memory functional in the limit of 
yt -C 1, first for one, and then for two fluctuators. Then 
the possibility of more than one jump is negligible, and 
we write for the probabilities of zero and one jump Pq = 
1 — yt and Pi = yt. The memory functional can be 
calculated by averaging over the time, ti, between the 
jumps, which at yt <C 1 are equally probable: 

P^ f* 7 

*i = P Q H / dh cos z/(t -ti) = 1- ~(vt - sinvt) . 

t Jo v 

If vt = (A/2Ejh)vt C 1 we can expand the sine and get 
SI/! = 1 — jis 2 t 3 /Q. At the optimal point v = and the 
fluctuator is not visible to the qubit. At other working 
points the dephasing rate is proportional to A 2 . 



Now take the case of two fluctuators. We have to cal- 
culate * 2 = (e^ 2 ^) where 



(12) 



Mt) = / d*>i(Xi-Xi) + ^(X2-X2) 
Jo 

+2A 12 ( X iX2-X?X§)] 
and A12 = V1V2/ E*h. Again, for y^t C 1 we get 



* a = i- 



±,i=l,2 



t _lj_ ( sin(^ ± \ 12 )t 



(13) 



Assuming v\, ^12 <C t 1 and expanding sines we get: 



71+72 ,2 , 3 



G 



n 1 



w±t 3 

6 



(14) 



We see that in this limit one can split the expression for 
^ in factors corresponding to the individual fluctuators 
just as in the case of linear coupling, but there appears 
an additional factor due to the nonlinear coupling. It 
is easy to see that a similar pattern will also appear for 
larger number of fluctuators. At longer times this prod- 
uct structure is lost. 

The interplay between the linear and quadratic cou- 
pling is now quite clear. If v A the linear coupling 
is dominant. Approaching the optimal point will reduce 
the dephasing until v becomes smaller than A where the 
term proportional to A 2 becomes most important. This 
contribution results from the interplay between the two 
fluctuators and cannot be eliminated; thus it represents 
the minimal dephasing possible at the optimal point. 
The physical reason for this is quite easy to understand. 
Switching of fluctuator 1 shifts the average point that 
fluctuator 2 is working around. Both positions of fluctu- 
ator 1 can not represent the optimal point with respect 
to fluctuator 2, and some dephasing is bound to occur. 
This is the most important physical insight that distin- 
guishes the quadratic coupling from the linear. In the 
case of quadratic coupling, even if the different fluctu- 
ators in themselves are independent, their effect on the 
qubit will be influenced by the positions of all the others. 
With linear coupling one finds that slow fluctuators, with 
yt -C 1 do not contribute to the dephasing. This is no 
longer true for quadratic coupling, as they play a role in 
determining the effect of the fast fluctuators even if they 
do not have time to switch during the experiment. Thus 
very slow fluctuators may be of great importance. 

b. Two-pulse echo: Perhaps more directly related to 
experiments are the echo signals. These are found using 
Eq. i|l(J|) where for two-pulse echo 



/?(*) = 



-1 for t < t , 

- 1 for r < t < 2t . 



(15) 



Here r is the delay between the initial pulse and the echo 
pulse, and the echo signal is centered around 2r. For 
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short times (jt < 1) this gives for one fluctuator 

2 7, 



\&i = 1 [yr — sinvr) 



and for two fiuctuators 



*2 = 1" E 

±,*=1,2 



27i^ - 7i 



sin(^ ± Ai 2 )t 

Vi ± Al 2 



(16) 



(17) 



IV. EXACT SOLUTION FOR SMALL NUMBER 
OF FLUCTUATORS 

We now turn to the calculation of the phase memory 
functional for arbitrary times. The method described 
here is in principle applicable to an arbitrary number 
of fiuctuators, but the resulting formulas quickly get 
impracticably large when the number of fiuctuators in- 
crease. The main idea is to consider the phase 4> as a ran- 
dom variable with some probability distribution p(cj>, t) 
that will depend on time. Once this is known the phase 
memory- function is given by 'J = J dx e^p(<p,t). By 
mapping to a correlated random walk problem we derive 
a Master equation for the probabilities p(4>, t) . The de- 
tails of how to calculate p(4>, t) are given in appendix lAl 



A. Distribution function for one fluctuator far 
from the optimal point 

Free-induction signal: Let us first discuss the results 
for the distribution function in the simplest case, that 
of one fluctuator far from the optimal point. The 
phase m emo ry functional for this problem was derived 
in Refs. 1 314 . and in the end we will rederive the same 



expression. However, it is the simplest example for il- 
lustrating the general method, and it gives new insight 
to the relation between the Gaussian approximation and 
the fluctuator model. To better understand the meaning 
of the results it is useful to recall the standard picture of 
dephasing by a Gaussian noise that is well known from 
NMR- physics (see, e. g., Ref . llOri . If the time t entering 
the phase memory functional 



$ = 



4>(t) = / dt'v(t') , 
Jo 



v(t) = v X (t) 



is much longer than the correlation time 
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of the 



fluctuating function v(t), the integral can be considered 
as the sum of a large number of uncorrelated contribu- 
tions. Consequently, by the central limit theorem, the 
phase will be distributed according to a Gaussian 



p(0) = 



1 



and the phase memory function is 



70 2 >/2 



From Eq. J^} for the correlation function we get 



(18) 



v 2 v 2 , 



2jt 



4-) 



t at jt > 1 . 



Thus, in the Gaussian approximation, the phase memory 
functional decays exponentially at jt 3> 1 with the rate 



.(G) 



v 2 /8j . 



(19) 



Using the method explained in Appendix lAl we can find 
an exact solution for the distribution function of <b 



J 



P(4>, t) = e 



- *-7* 



. T it/ 2 ) + - 
v 



{t±24>/v)h(i^/t 2 ~(2^/vY 
I 



v 



2 A _ t 
v 



(20) 



where the different signs correspond to different initial 
states of the fluctuator and I\{z) is the modified Bessel 
function. Without jumps of the fluctuator, the result 
would be only the moving <5-pulse of a constant ampli- 
tude, 6 (cf> =F vt/2). The value vt/2 is the maximal possi- 
ble value of 4> acquired for the time t, while the jumps of 
the fiuctuators account for the smooth part. Averaging 
over the initial state of the fluctuator we get (not writing 



the 6-1 unctions) 



P(<M) 



-7* 



vt/2) + 8{4> + vt/2) 



7_I 1 ( 7 i v /l 
v 



(2/vt) 2 4> 2 ) 



Vi - fflWP 



(21) 



This is plotted in Fig. [3 for the times t = 1, 5, 10, 7 = 1 
and v = 1. We observe that the central region is similar 
to a Gaussian, but at short times this is cut off by the 
(5-functions represented by the vertical lines. At yt ^> 1 
the function is indeed close to a Gaussian, as can be seen 
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8 10 
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FIG. 2: Distribution function at time 1, 5, 10 for 7 = 1 and 
v = 1. Inset: time dependence of the integrated strength of 
the central peak. 



from the asymptotic expansion of the Bessel function 



-it 



1 



j 1 ( 7 yi-(2HV) ^ 

2y/l - (2/i/i) V ~ v^T*' 



Comparing to Eq. I|18|l we see that we recover the result 
for the dephasing rate in Eq. lj*H3|). 

However, we know from Ref. |^ that if v > 27 we have 
pronounced non-Gaussian behavior. We can now under- 
stand this from the point of view of the distribution func- 
tion. The smooth central part of this indeed approaches 
a Gaussian for jt >■ 1 and this gives the dephasing rate 
(|19|l . but the 5-functions at the ends only decay at the 
rate 7. As long as the r^ G ' of l|19fl is smaller than 7 
the decay will be controlled by the central part and the 

(G) 

Gaussian approximation is valid. If 1^ > 7 the decay 
is limited by the (5-functions, and is set by the rate 7. 

Using the distribution function 1|21[1 one can calculate 
the phase memory functional (see Appendix [SJ and one 
finds for 7< <C 1 an exponential decay with rate 



r, 



7 -Re (V7 2 ~ f 2 / 4 ) 



(22) 



= 1 for the 



Figure [3| shows as function of v at 7 
Gaussian and the nuctuator models. 

Two-pulse echo: The same method as outlined in the 
Appendix [S] allows one to calculate the phase distribu- 
tion function, p e (0, r) for the two-pulse echo signal. The 
result can be expressed in the form 



Pe 



7 / „ cos fc^siniUfcT (7 sinwfcT + Wk cos Wkr) 
-— / dk ^ 



Here Wk = \J (kv/2) 2 — j 2 , while r is the delay time 
between the first and second pulse. The smooth part of 
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FIG. 3: Dephasing rate as function of v, 7 = 1. The solid 
line is the exact result for one nuctuator, Eq. (1221 . and the 
dashed line is the Gaussian approximation, Eq. 1191 . 



the distribution given by the second item in the above 
formula is plotted in Fig. It is qualitatively similar to 
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FIG. 4: Distribution of the phases for the echo signal for v = 7 
and different products 7T (shown near the curves). 

the phase distribution for the free induction. 



B. Dephasing rate as function of working point 

Let us then consider the decay of the phase memory 
functional for different working points . Again we recall 
the situation in NMR-physics where the loss of the signal 
after the spins are set precessing by a ir/2- pulse is caused 
by two independent processes. The phase memory func- 
tional considered above measures the random contribu- 
tions to the phase caused by fluctuating energy difference 
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between the two states (which in an NMR experiment is 
caused by fluctuations in the magnetic field parallel to 
the external main field). In addition there are processes 
which flip the spin from one state to another, so called 
Xi-processes. These are caused by fluctuations in the 
magnetic field normal to the external field. If we denote 
the decay rate of the excited state into the ground state 
Ti and add the two contributions we have the total decay 
of the spin precession signal 



1 



(23) 



The factor 1/2 in front of Ti can be understood from 
the fact that if the probability of the excited state de- 
cays with rate r 1; the amplitude decays with the rate 
Fi/2 and this is what enters the off-diagonal elements of 
the density matrix. For a more detailed discussion see 
Ref. ITol or, for qubits, Ref. 

Let us now discuss how the relative strength of the 
two terms of Eq. I|23|) changes as we change the working 
point of the qubit. Looking back at the Hamiltonian @ 
we remember that it is equivalent to a spin 1/2-particle 
in a static magnetic field B = Eje x + Ae z while the 
noise is always along the z-axis, v — ve z . We denote the 
angle between B and the z-axis by 6 = arctan(_Ej/A). 
In particular, 9 = corresponds to working far from the 
degeneracy point where <5 3> Ej while 9 = tt/2 is the 
degeneracy (optimal) point A = 0. The time evolution 
of the qubit is then a precession on the Bloch sphere 
around the total field (Fig.[SJ). Far from the degeneracy 



Noise field 



Linear regime 
Main field 




FIG. 5: 

sphere 



Degeneracy point 



Different working points represented on the Bloch 



point the main field is directed along the z-axis (pointing 
at the north pole of the sphere), whereas at degeneracy 
it is along the £-axis (on the equator). All the time the 
noise vector v is in the z-direction. We see that the noise 
component parallel to the external field, which gives the 
r^, is largest far from degeneracy, which agrees with our 
previous discussion. However, the noise normal to the 
field, giving Ti-processes, is maximal at degeneracy. In 



the Gaussian approximation this is given by (see Ref. E3) 



dr (v(t)v(t + t)) cos 



Using the correlator JSJ one can rewrite this expression 
for the case of N identical fluctuators with parameters v 
and 7 as 



2 El 



4 7 2 



El 

Eo 



(24) 



The Gaussian approximation for the is more difficult 
to obtain because the square root in the energy J7J) makes 
the average in Eq. (|18|l not treatable analytically. How- 
ever we can expand in the lowest order in v/Eq, which is 
a good approximation except a close vicinity of the de- 
generacy point, where the coefficient in front goes to zero 
and higher order terms need to be calculated. This gives 
the same result as the rf d = 1/T 2 ad of Ref. OH 



.(G) 



-cos 2 6» / dT{u(t)u(t + T)) , 



(I 



which for N identical uncorrelated fluctuators yields 

2 



(G) ^Nf_ /A. 
1 7 V£o 



(25) 



Now we want to compare these Gaussian results with the 
exact expressions found using the method of Appendix 
lAl In FigureOthe relaxation times (inverses of the relax- 



In T 




FIG. 6: Relaxation times for the case of two fluctuators as 
function of working point for weak coupling, v/Ej = 0.1, and 
j/Ej = 0.1. A is measured in units of Ej. 1 - 1/r^ , 
Eq. 2-Ti, Eq. JHJ, 3 - T 2 = (l/2Ti + 1/T^)" 1 , 4 - T$, 
calculated according to exact expressions from Appendix 1X1 
5, 6 - results of numerical simulation for T\ and T2, respec- 
tively 

( G) 

ation rates) T\ (curve 2) and 1/rV (curve 1) are shown. 
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They are calculated according to Eqs. {^ty and lj5Sfr. re- 
spectively. The decay time, — ri 1 (curve 4), of the 
phase memory functional calculated using the method of 
Appendix lAland the resulting decay time, T 2 (curve 3), of 
the spin signal according to Eq. I|23|l are also shown. The 
points represent the rates obtained by a numerical simu- 
lation of the time evolution according to the Hamiltonian 
(|2*|). which performs averaging over many realizations of 
the random process. All curves are calculated for the case 
of 2 fluctuators with coupling strength v/Ej = 0.1 and 
switching rate 7/Ej = 0.1. Thus v/7 = 1 and this case 
belongs to the so-called weak coupling regime^ We see 
that for all working points the decay is well described by 
the Gaussian approximation. This is because the rates 
always are slower than the limiting rate 7 set by the cor- 
relation time of the fluctuators, similar to what was de- 
scribed above for the case far from the degeneracy point. 
For most working points the rate 3> Ti and this dom- 
inates the T 2 . Close to the degeneracy point we see that 
Ti becomes more important and at degeneracy it domi- 
nates completely giving T 2 ~ 1^/2. The same quantities 



InT 




0.2 0.4 0.6 0.8 1 



A 

FIG. 7: Relaxation times as function of working point for the 
case of strong coupling. Two fluctuators with v/Ej — 0.1, 
'y/Ej = 0.01. The legend is the same as in Fig. |S| 

for the case of strong coupling v/Ej = 0.1, 7/E'j = 0.01, 
u/7= 10 are shown in Figure^] Far from the degeneracy 
the situation is similar to the one described earlier, with 
the rate, 7, determined by the decay of the (5-function 
peaks in the distribution function. Closer to the degen- 
eracy the rate is slower and the Gaussian approximation 
gives good results. 

Note that in the case of strong coupling the difference 

(G) 

between the approximate expression and the exact 
becomes more noticeable. At large A this is because 
of the essentially non-Gaussian character of the noise, as 
discussed before. At small A the difference comes from 

( G) 

the fact that the only gives the Gaussian approxima- 
tion to r^, in the lowest order in v/Ej. We expect that 



if we had calculated the phase memory function in the 
Gaussian approximation according to (|18|) and using the 
exact expectation value ((f> 2 ) (which is hard to find ana- 
lytically) the result would agree completely with the 
calculated by the method of Appendix lAl since the decay 
time is much longer than the correlation time 1 /j of the 
noise and the central limit theorem should work. Note 
also that at degeneracy the decay is still dominated by 
the phase relaxation processes, T^, while the Ti-processes 
only give a small correction. 

Figures El and should be compared to the experimen- 
tal results of Astafiev et almL There is clear qualitative 
agreement, but it is hard to try to make a quantitative 
fit, especially for the T 2 where there is very little data. 

Let us rather look back at the Hamiltonian J5J) and Fig- 
ure |3J So far we only considered noise in the a z compo- 
nent of the Hamiltonian, as appropriate for noise sources 
coupled to the charge of the qubit. As discussed in Ref.0 
there is also the possibility of noise in the Joscphson cou- 
pling (the a x part). In Figure [3] this would correspond to 
a noise vector parallel to the x-axis. This noise would be 
transversal far from degeneracy, giving large T 2 and small 
T\ and it would be longitudinal at degeneracy point, with 
Ti large and T 2 small. The fact that the experimental 
results are similar to our figures HH and [7| rather that the 
opposite shows that in these experiments the noise in the 
<j z part is dominant. 



V. DISCUSSION 

At the end we would like to discuss two simple obser- 
vations based on the above results. 



A. Measuring the distribution function? 

The method that we have used to find the phase mem- 
ory functional gave as an intermediate result an expres- 
sion for the distribution function of the phase p((f>, t). The 
phase memory functional is the average of the quantity 
e 1 ^ with this distribution. But can the full distribution 
function be compared to experiment? In the first set of 
experiments^ only averages could be measured, but re- 
cently there has been demonstrated single shot readout 
of the qubit stated Is it possible to extract the full dis- 
tribution function from such data? Unfortunately, the 
answer is no for the following reason. For each realiza- 
tion of the experiment, there is a certain realization of 
the random external noise. This gives the phase differ- 
ence 4> a certain well defined value. The next realization 
gives another realization of the noise process and another 
value for the phase difference and so on. Since we have 
no knowledge about which realization of the noise is rel- 
evant to a certain experiment, we have to describe the 
final state of the qubit by a density matrix representing 
a mixed state. Since all outcomes of all possible exper- 
iments can be calculated from the density matrix, no 
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more information about the system can be learned that 
what is contained in the density matrix. Since the phase 
distribution was calculated ignoring the relaxation (Ti) 
processes the density matrix will be represented by a vec- 
tor in the equatorial plane of the Bloch sphere (but for a 
mixed state it will not be on the surface of the sphere but 
at some interior point), and the components of this vec- 
tor are exactly the real and imaginary parts of the phase 
memory functional. Thus, the phase memory functional 
contains all information that we can extract about the 
qubit through experiments. A different way to express 
the same is that there are many different phase distri- 
bution functions yielding the same qubit density matrix, 
and there is no way that one can experimentally distin- 
guish these in an experiment with a single qubit. 



B. Non-Gaussian Ti? 

Notice that it seems from the figures El and that the 
energy relaxation time T\ is always well described by the 
Gaussian approximation. Can we understand this in a 
better way? The Gaussian approximation is good pro- 
vided the correlation time of the noise is much shorter 
than the decay time, T± or T2. In that case separation 
of timescales enables one to deduce an exponential decay 
with ratesti^ 

l/Ti oc S(E ), 1/T 2 = l/2Ti + Tf, Tf oc S(0) 

where S(u) is the noise power spectrum. For the 
telegraph process with switching rate 7 and coupling 
strength v we find that the time correlation is given by 



(xWx(o)) 



_ p -27l*l 



so the correlation time is I/27 and the noise power spec- 
trum 



This gives the relative rates 



7 2 + uj 2 



1 



T17 7 2 + 



< 



v 



1 2 

(X 

7 V7 



We see that as long as the noise is weak compared to the 
qubit splitting, v < Eq, the condition 7T1 3> 1 for the va- 
lidity of the Gaussian approximation, is always satisfied. 
For the T2 the situation is different, and the condition 
T 2 < 7 is violated if v > 7, leading to non-Gaussian be- 
havior—. So in the case of dephasing the Gaussian approx- 
imation so to speak predicts its own breakdown whereas 
for energy relaxation it is consistent as long as the noise 
is weak compared to the level spacing. 
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APPENDIX A: CALCULATION OF 
DISTRIBUTION FUNCTION 

1. Single fluctuator 

To explain the method we re-derive the result for one 
fluctuator with linear coupling since this is the simplest 
case. Let us for simplicity also assume that f3(t) = 1 at 
t > and start with the calculation of ^i(t). 

Let us discretize the integral (|T2")) for A i2 = introduc- 
ing small time steps r = t/N, where N 3> 1. Then the 
random phase shift cf>(t) can be expressed as 



N 



VT 



Xr. 



n=l 



Hence, the integration over time can then be thought of 
as a random walk process, where at each time step the 
random walker moves a step a — tv/2 in the direction 
depending on the current position of the fluctuator. The 
steps are correlated, but only with the previous step. The 
probability that a step is in the same direction as the 
previous one is a = 1 — jt and the probability for a step 
to be in the opposite direction is (3 = jt. Let m be 
the number of steps from the origin (so that the position 
is x = am). We want to find the probability P n (m) 
to be in position m at time step n (dimensional time 
t n = nr). This is found by the following method. We 
split the probability in two parts: the probability to reach 
point m coming from the right, A n (m) and from the left 
B n (m), so that P n (m) = A n (m) + B n (m). We then have 
the equations 



A n+ i(m) = aA n (m - 1) + f3B n (m - 1) 
B n +i(m) = f3A n (m + 1) + aB n (m + 1) 



(Al) 



We need the continuum limit, letting N — * 00, and t — > 
with Nt = t fixed. Writing 



a(4>, t) = a(m<j, tit) = A n (m) , 
b(4>,t) = b(mcr,nT) — B n (m) 

and expanding to first order in r we get 

a + rat = a{a — aa^) + j3(b — ab^) , 
b + rb t — (3(a + aa^) + a(b + ab^) . 



(A2) 



(A3) 



Here subscripts <f> and t denote partial derivatives with 
respect to <f> and t, respectively. Adding and subtracting 
these we get (with p = a + b and q = a — b) 



rp t = {fi- a)aq<t, , 
q + rq t = (a - /3)q - o-p^ . 



(A4) 



Differentiating the second of these and inserting q^ from 
the first we obtain the final equation for p 



Ptt 



V<t>4> 



(A5) 
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which is called the telegraph equation. We guess the 
solution p — e J ( fc *-^*) arK } get the dispersion relation 



~2 ry2 



(A6) 



w± = -«7±y(^) « 2 

The general solution is then 

^ [a K e--+* + 6 K e---*] e te * . (A7) 

-oo ^ 7r 

The coefficients {a K , 6 K } can be obtained from the initial 
conditions, 

/oo 
dn [a K + b K ] 
-CO 

a K + b K = 1 , (A8) 

/OO 
d/c [o; + a K + u)-b K ] e lkx 
-co 

==>• uj + a K +u-b K = -KVXo , (A9) 



which yield 



2 2^kV 2 /4- 7 2 



(A10) 



This is to be inserted into Eq. i|A7(l . The result is given 
by Eq. (JUJl. 

The memory functional can be expressed as the expec- 
tation value *i = / d(j)p{4>, t) e^. Thus 



2[i 



dn [a K e-^ +t + 6«e _<w -*] <5(k + 1) 

f 



(All) 



where /i = y/T— v 2 j^i 2 - This agrees with the result of 
Ref . H 



Two fluctuators 



Now we turn to the case of two fluctuators. Again we 
discretize Ijl2(l to get 



N 



There are now three kinds of steps, depending on the 
settings of the fluctuators. If both have the value +A 
there is a step a = v + -|, if both are — \ there is a step 



and for one of each the step is 7 



We 



have the following probabilities for each jump, depending 
on the previous state: 



Previous ■ 
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■ 2 2a(3 (3 2 
7 a/3 a 2 + (3 2 a/3 

(3 2 2af3 a 2 
The total probability must now be split in 3 parts 
P n {4>) = A n {4>) + B n {4>) + C n (<f>), where A n {4>) is the 
probability to reach point (j) a t time step n with a a 
jump, B n (4>) with a (3 jump and C n {4>) with a 7 jump. 
We have then the set of equations 



A n +i(4>) = a 2 A n (ct) - or) + 2 B n (4> - dr) + a/3C n ((f> - dr) , 
B n+1 (cf>) = p 2 A n {cp - p T ) + a 2 B n {4> - pr) + apC n [cj> - pr] , 
C n+1 (4>) = 2aPA n (0 - 7 r) + 2a$B n {cp - 7 r) + (a 2 + p 2 )C n ( 

I 



- 77-) . 



Again we introduce continuous variables a, 6, c and ex- 
pand to first order in r to get 



a t — —27a — aa0 + 7c , 
b t = -276 - fibj, + 7c , 
ct = 27a + 276 — 27c — 7C0 . 

This can be written in matrix form 



a t = Ma, 



(A13) 



(A14) 



(A12) 



where 

(-27 - ddj, 7 

-2 7 -/39 7 I . (A15) 

27 2 7 -2 7 + 7 <9 

We then guess the solution in the form 

a = Ae*^ - ^ (A16) 

which gives the eigenvalue equation 

-iwa=MA (A17) 
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where M is the matrix M with 9^ replaced by in. From 
this we get the dispersion equation for u> in terms of k: 



to 



6ji — k— ] jj 1 



1 2 
- K V 



87 2 - Aqin^ - K 2 y ) x 



-■2- i;r l> 2 + y 



3 A 
K — I V 

2 



2 -l=o. 



This equation has the three solutions: ujq and u)± where 
ujq is the solution that goes continuously to — 2ij when 
A — > 0. The general solution is then 

a = / p A K e -*M , (A18) 

where 



T K 



(A19) 







( a K ^ 






is a matrix of coefficients that has to be determined by 


a K = 








the initial conditions. 











functions a(^>, 0)... and the first two derivatives. How- 
ever, because of the special form of the equations we 
can calculate all derivatives at t — from the functions 
a(0,O)... 

at(0,O) = Ma(^,0), a tt (0,O) = M 2 a(0,O) . (A20) 

The typical initial conditions would then correspond to 
specifying the initial type of jump in the random walk. 
For example if this was of type A we would have a(x, 0) = 
S(x) and b(x, 0) = c(x,0) = 0. Note that in this simple 
case where only one of the a, b and c are nonzero at t = 
the procedure could be simplified by writing the general 
solution for p = a + b + c and initial conditions for this. 
The more general case would be that all of a, b and c 
are nonzero and the complete matrix A K is needed. This 
would be the case for example when calculating echo sig- 
nals, where the equations after the echo pulse has to be 
solved with initial conditions of this type corresponding 
to the solution of the equations before the echo pulse is 
applied. 

Introducing the Fourier transformed functions (in <ft) 



(A21) 



Since we have three coefficients to determine for each of 
the a, b and c it appears that we need to specify both the 



we can write the initial conditions as 



J 



a*(f = 0) =A K | 1 



Ma K (t = 0) = -iA K 




M 2 a K (t = 0) = -A K 



(A22) 



These can be written more compactly if we introduce the matrix M a with the left hand sides of the above equations 
as columns 



1 -iw (~iu} ) 2 

( a K (t = 0) | Ma K (t = 0) | M 2 a K (t = 0) J = A K | 1 H^+) 2 

1 — icij_ (— icj_) 2 



(A23) 



m„ 



r 



From which the coefficients are found as A K = M a Q . 
The final solution is then 



P(<t>,t) 



' d, K r 



iujQt 



+ b K e- iu + f + 



c K e e v 



where a K — a J. and similarly for b K and c K . Again the 
average is calculated from 



<e*> 



dt/M&pfa t) 



a _ ie -i^ot + b _ ie -iu+t + c _ ie - 



where k = — 1 in the u> in the last expression because of 
the 5-function from the <f> integral. 

Let us find the explicit expressions for the coefficients 
a K , b K and c K in this case. Adding the lines in Eq. (|A23|) 
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we get 

~\~ b K ~\~ c K 1 j 
a K uj + b K u)+ + c K U- = k(v(xi + xl) + 2^X1X2) = A K 
a K wl + b K uj\ + c K u 2 - = « 2 (Kx? + xl) + 2Ax?X°) 2 

-2i 7 /«/(x° + Xl) - 8»7«AxSx2 = B * ■ 

Here x? 2 represent the initial state of the fluctuators. 
Also of interest are the values of these averaged over the 
initial states of the fluctuators. Assuming all 4 settings 
arc equally probable we have 



0, B° 



In terms of these the coefficients are expressed as 
(w - A){ujq + uj+) - (ujI - B) 



b K 



(A24) 



a K = 



(co - U>-)(U)+ - UI-) 

-(luq - A)(lo + uj_) + {ul - B) 

(u) - uj + ){u + - w_) 
1 - b K - c K . 

3. General case 



The above method is in principle simple to generalize 
to any number of fluctuators, but the number of equa- 



J 



tions increases exponentially in the number of fluctua- 
tors. 



The general equation is 



a t = Ma, 



a = 



b 

c 

V : J 



Guessing the solution a = Ae^ K ^~"^ we get the dis- 
persion equation — iuja — MA, which determines the n 
eigenvalues LOi(n) (i = l...n) as functions of k. Here 
n = 2 N with N the number of fluctuators. The general 
solution is 



/ 



2^ 



( 



-iu>2t 



\ 



Defining 



a K — 



/ a K 
K 



A K e u , and M = M{d$ -> in) 



we determine the coefficient matrix A K from 



(a K (i = 0)|Ma K (t = 0)|--- | M n - x k K {t = 0) ) = A K 



M a 



( 1 -iW! • • • (-^i)™" 1 \ 



\ 1 -iu n ■ ■ ■ (-iu n ) n 1 / 



The matrix M K can be written as the contraction of a 
third order tensor Mt with a K . Writing the tensor indices 
we have 

[M„]ij = [M T ]«fc[a«(0)]fc, [M T ]ijk = 
We then get 

[A K ]ij = [AfT] 4 mfcP _1 ]mj[a K (0)] fc = [M Q ] ijk [a K (0)] k 
and 

M*)]i - [AJyrt - [T«]ife[a K (0)] fe , 
[T K ] ik = [Mnlo-fcrt-. 

One can also write equations for echo experiments. Con- 
sider the situation where we initially prepare a state, then 
apply the echo pulse at time t e and then measure the 
state at the time 2t e when the echo signal appears (two 



pulse echo). The state just before the application of the 
echo pulse has to be calculated as above and then this is 
used as the initial state for the evolution after the echo (it 
is assumed that the duration of the echo pulse is short). 
After the echo pulse the matrix M is changed because all 
jumps of the random walk changes sign. Then also the 
uji change. Let M~ , iv~ represent quantities before the 
echo pulse and M + , uof after the pulse. Then 

[a«(te)]i = [T-}ik[a K (0)] k 
[a K (2t e )]i = [T+] ik [a K (t e )] k 

with 

[T~] ik = [M^ ijk [e^}j, [T+] ik = [M+] ijk [e u+t ']j 
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